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ABSTRACT 

A multisite photometric campaign for the /3 Cephei and eclipsing variable 16 Lacertae 
is reported. 749 h of high-quality differential photoelectric Stromgren, Johnson and 
Geneva time-series photometry were obtained with ten telescopes during 185 nights. 
After removing the pulsation contribution, an attempt was made to solve the resulting 
eclipse light-curve by means of the computer program EBOP. Although a unique 
solution was not obtained, the range of solutions could be constrained by comparing 
computed positions of the secondary component in the Hertzsprung-Russell diagram 
with evolutionary tracks. 

For three high-amplitude pulsation modes, the uvy and the Geneva UBG ampli¬ 
tude ratios are derived and compared with the theoretical ones for spherical-harmonic 
degrees € < 4. The highest degree, £ = 4, is shown to be incompatible with the ob¬ 
servations. One mode is found to be radial, one is £ = 1, while in the remaining case 
€ = 2 or 3. 

The present multisite observations are combined with the archival photometry in 
order to investigate the long-term variation of the amplitudes and phases of the three 
high-amplitude pulsation modes. The radial mode shows a non-sinusoidal variation on 
a time-scale of 73 yr. The ^ = \ mode is a triplet with unequal frequency spacing, 
giving rise to two beat-periods, 720.7 d and 29.1 yr. The amplitude and phase of the 
£ = 2 or 3 mode vary on time-scales of 380.5 d and 43 yr. 

The light variation of 2 And, one of the comparison stars, is discussed in the 
Appendix. 

Key words: stars: early-type - stars: individual: 16 (EN) Lacertae - stars: individual: 
2 Andromedae - stars: eclipsing - stars: oscillations 
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1 INTRODUCTION 

16 (EN) Lacertae = HR 8725 (B2IV, V = 5.59), a mem¬ 
ber of_I^c^DBlaj_Js_a^_single;lmed_spectroscopic binary 
dLed I 19 IOI : IStmve fc Bobrovnikofil Il92^ and an eclipsing 
variable. The orbital period, derived from the epochs of 
minimum light, Pprb = 12.09684 d Jjerzvkiewicd Il980l : 
IPigulski fc Jerzvkiewicj Il988l l. The system consists of a 
well-known P Cephei variable and an undetected secondary. 
The P Cephei variation of the primary is dominated by three 
pu lsation mode s with frequencies close to 6 d“^. According 
to iFitchl ll 19691 . henceforth F69), who based his analysis on 
archival 1951, 1952 and 1954 radial-velocity and photoelec¬ 
tric blue-light observations, the first two modes have con¬ 
stant amplitudes, while the amplitude of the third mode 
varies on a time scale of months. F69 described the hrst 
two modes as singlets, having frequencies /i = 5.91134 and 
/2 = 5.85286 d“^, and the third mode as a doublet, con¬ 
sisting of two terms with frequencies /sp = 5.49990 and 
/a ,2 = fi — 5/orb = 5.49799 d“^. Comparable frequencies 
were derived from the 1964 and 1965 photoelectric obser¬ 
vations, obtained by one of us (MJ) at Lowell Observatory, 
viz. /i = 5.91120 ± 0.00005, f 2 = 5.85503 ±0 00010 , and 
/a = 5.50322 ±0.00009 d'^ lljarzebowski et al.lll979l'l . The 
hrst value agrees with that of F69, but the other two dif¬ 
fer from F69’s /2 and / 3.1 by 0. 0022 and 0.0033 d~\ re - 
spectively. From these differences, Ijarzebowski et ^ ill979l l 
concluded that F69’s values of /2 and /a,i suffered from an 
error of one cycle per year (yr“^). 

An analysis of all photometric observations of 16 

Lac obtained through out _ 199 2 wa s carried out by 

Ijerzvkiewicz fc Pigulskil ll 19961 . Il999l . henceforth JP96, 
JP99). The main results of the analysis can be summarized 
as follows: ( 1 ) the amplitudes of the /i and /2 modes vary on 
a time-scale of decades, the reciprocal time-scales amount¬ 
ing to 0.014 and 0.020 yr~^, respectively, (2) the third mode 
is conhrmed to be a doublet, but with frequencies different 
from those derived by F69, viz. /s^ = 5.5025779 ± 0.0000005 
and /a ,2 = 5.5040531 ±0.0000008 d“^. Note that none of 
these frequencies bears a simple numerical relation to the 
orbital period. Moreover, JP96 demonstrated that there is 
no correlation of the pulsation amplitudes with the orbital 
phase. 

The radial-velocity (RV) data available at the time were 
shown by JP96 to be consistent with the above-mentioned 
photometric results. The RV data, however, were much less 
numerous than the photometric data, making this conclu¬ 
sion somewhnt_UMertato. The situation has improved af¬ 
ter iLehmann et al.l (1200 ll . henceforth LOl), provided new 
RV observations, more than doubling the number of exist¬ 
ing measurements. A periodogram analysis of these and the 
older data led LOl to a number of frequency solutions, the 
details of which depended on the weights assigned to the 
RV measurements from a particular source. As far as the 
periods are concerned, the results valid for all weighting 
schemes can be summarized as follows: ( 1 ) the hrst mode 
is an equidistant triplet with the central term having the 
largest amplitude and a period. Pi, equal to 0.16916707 d 
(/i = 5.9113160 d“^). The remaining two periods, 
and Pp in the notation of LOl, are equal to 0.16916605 
and 0.16916809 d, respectively. The frequency separation of 
the triplet amounts to 0.0000356 d”*^, corresponding to the 


amplitude-modulation frequency of 0.0130 yr”*^, (2) the sec¬ 
ond mode is a doublet; in the order of decreasing amplitude, 
the periods, P 2 and P 2 in the notation of LOl, are equal to 
0.17085553 and 0.17077074 d, respectively (/z = 5.8528981 
and = 5.8558041 d“^), resulting in a beat period of 
344 d (the beat frequency of 0.00291 d~^), (3) the third 
mode is also a doublet; in the order of decreasing ampli¬ 
tude the periods are P 3 = 0.18173251 and P^ = 0.18168352 
d (/s = 5.5025928 and = 5.5040765 d“^); in this case 
the beat period is equal to 674 d. Conclusions (1) and (3) 
approximately agree with the results of JP96 and JP99, 
but conclusion (2) does not. The disagreement is twofold. 
First, the time-scale of the amplitude modulation of the 
/2 term derived by LOl is less than one Hftieth of that 
derived by JP96 and JP99. Second, LOl’s / 2 , i.e. the fre¬ 
quency of the higher-amplitude component of the / 2 , 
doublet has a value c lose to that originally de termined by 
F69 and dismissed by Ijarzebowski et al.l (Il979h as a yearly 
alias of the photometri c (2 v alue they derived. Note that 
/2 of ijarzebowski et al.l lll979ll is close to f^, the frequency 
of the smaller-amplitude component of the / 2 , P 2 doublet. 
This, of course, is the consequence of the doublet’s beat- 
period having its value close to 1 yr. 

In addition to the three modes just dis cussed, six fainter 
terms were detected by I Jerzvkiewicj (Il993fl in his 1965 data. 
The y amplitudes of these terms ranged from 2.1 ±0.14 to 
0.7 ±0.14 mmag, and the frequencies were equal to (in the 
order of decreasing amplitude) 0.1653, 7.194, 11.822, 11.358, 
11.414, and 11.766 d~^. The first of these is equal to twice 
the orbital frequency, suggesting ellipsoidal variability. How¬ 
ever, the observed amplitude and phase excluded this possi¬ 
bility. The third frequency is equal to 2/i and the three last 
frequencies are the combination terms /2 ± fa . fi ± fs and 
fi ± / 2 . The 7.194-d~^ term was attributed by I Jerzvkiewicj 
h993ll to an independent pulsation mode and was used 
as such, together with the th ree strongest ones, in an as- 
teros eismic study of the star (iDziembowski fc Jerzvkiewicj 
199611. Subsequently , how ever, this term was shown by 


Sarevan et al.l (119971 1 and iHandler et al.l (l2006fl to be due 
to a li ght variation of 2 Andromedae, used by Ijerzvkiewic j 
1I1993II as a comparison star. 

During a multisite photometric campaign carried out 
between 2 August 2003 and 9 January 2004, 16 Lac has been 
observed together with the P Cephei variable 12 (DD) Lac¬ 
ertae. Results of the observations of 12 Lac and their an aly- 
sis were published some years ago llHandler et al.|[2006ll . In 
the next section, we describe the 2003-2004 multisite cam¬ 
paign’s photometric observations and reductions. In Section 
3, we carry out a frequency analysis of the campaign’s uvy 
time-series of 16 Lac. In Section 4 we use the analysis re¬ 
sults to remove the intrinsic component of the variation from 
the time-series, thus bringing out the eclipse, and discuss the 
eclipse solutions and the evolutionary state of the secondary 
component. In Section 5, we derive the primary component’s 
fundamental parameters. Section 6 is devoted to determin¬ 
ing the harmonic degree of the three highest-amplitude pul¬ 
sation modes of 16 Lac. The long-term variation of the pho¬ 
tometric amplitudes and phases of these three modes is in¬ 
vestigated in Section 7 using the present and archival data. 
The last section contains a summary and discussion of the 
results. The light variation of 2 And is examined in the Ap¬ 
pendix. 
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2 OBSERVATIONS AND REDUCTIONS 


Our 2003-4 photometric observations of 16 Lac were carried 
out at ten observatories on three continents with small- to 
medium-sized telescopes (see Table 1). In most cases, single¬ 
channel differential photoelectric photometry was acquired. 
At the Sobaeksan Optical Astronomy Observatory (SOAO) 
and the Bialkow Observatory (BO) the photometry was 
done with CCD cameras. Wherever possible, Stromgren uvy 
filters were used. At the Sierra Nevada and San Pedro Martir 
(SPM) Observatories simultaneous uvby photometers were 
available, including b filters as well. However, the u data from 
SPM were unusable. At BO, a Stromgren y filter was used. 
At four observatories where no Stromgren filters were avail¬ 
able we used Johnson’s V. Finally, as the photometer at the 
Mercator telescope had Geneva filters installed permanently, 
we used this filter system. We chose the two ‘classical’ com¬ 
parison stars: 10 Lac (09V, V = 4.88) and 2 And (A3Vn, 
V = 5.09). A check star, HR8708 (A3Vm-|-F6V, V = 5.81), 
was additionally observed during one of the SPM runs. 
HD 216854 (F5, V = 7.31) was used as a comparison star at 
SOAO; its constancy was checked against two fainter stars. 
At BO, BD-b40°4950a = PPM 63607 = CSC 3223-01835 
(F5/K0, V = 9.29) was used as a sole comparison star. In 
order to compensate for the large brightness difference be¬ 
tween the program and the comparison star (~3.7 mag), BO 
observers took a sequence of CCD frames with short and 
long exposure-times, such that 16 Lac was well exposed on 
the frames with short exposure-times, while the comparison 
star, on the frames with long exposure-times. Reductions 
of the SOAO and BO observations included calibrating all 
CCD frames in a sta ndard way and processing them with the 
Daophot H package dStetsorJ 19871 '). Then, magnitudes were 
obtained by means of aperture photometry, the magnitudes 
of the comparison stars were interpolated for the times of 
observation of 16 Lac and differential magnitudes were cal¬ 
culated. In the case of the BO data, several consecutive data 
points were averaged, resulting in the final time-series. 


As we mentioned in the Introduction, 2 And turned out 
to be a low-amplitude periodic variable. We shall discuss 
the variability of 2 And in the Appendix. No evidence for 
photometric variability of 10 Lac was found. We thus pro¬ 
ceeded by pre-whitening the variability of 2 And with a fit 
determined from all its differential magnitudes relative to 
10 Lac. The residual magnitudes of 2 And were then com¬ 
bined with the 10 Lac data into a curve that was assumed 
to reflect the effects of transparency and detector sensitivity 
changes only. These combined time-series were binned into 
intervals that would allow good compensation for the above- 
mentioned non-intrinsic variations in the target star time- 
series and were subtracted from the measurements of 12 and 
16 Lac. Note that the binning minimizes the noise in the dif¬ 
ferential light curves of the targets. Finally, the photometric 
zero-points of the different instruments were compared and 
adjusted if required. In particular, adjustments were neces¬ 
sary for the SOAO and BO observations because they were 
obtained with different comparison stars. Measurements in 
the Stromgren y and the Johnson and Geneva V filters were 
analysed together because these filters have very nearly the 
same effective wave-length; the combined y, V light-curve is 
henceforth referred to as ‘the y data.’ For further details of 


the re ductions, common to 12 and 16 Lac, see [Handler et aP 

(l2006tl . 

For both stars, 12 and 16 Lac, the y data were the most 
extensive by far. In the case of 16 Lac, there were 3190 y- 
Hlter measurements (including Johnson and Geneva U-filter 
measurements, see above), 2012 u-filter measurements, and 
1686 w-filter measurements. The Mercator, Fairborn, Lowell, 
and Mayaki data included measurements falling between the 
hrst and the fourth contacts of six eclipses. Omitting these 
measurements resulted in the time-series suitable for fre¬ 
quency analysis; in the process, the number of measurements 
was reduced to 3055, 1895, and 1583 in y, v, and u, respec¬ 
tively. In all cases, the measurements spanned an interval of 
179.2 d. The number of measurements in the Geneva filters 
was 411, with 395 taken outside eclipses. In both cases the 
data spanned 129.9 d. A sample of the light curves of 16 Lac 
is presented in Fig. [T] 


3 FREQUENCY ANALYSIS 

Using the results of JP99’s analysis (see the Introduction) 
we can predict blue-light amplitudes of the fi {i = 1, 2, 3) 
modes for the epochs of our 2003-4 observations. In the case 
of /i and / 2 , the amplitudes predicted under an assumption 
of sinusoidal variations with the reciprocal time-scales of 
0.014 and 0.020 yr”*^ amount to 12.4 and 15.1 mmag, respec¬ 
tively, and both amplitudes should be very nearly constant 
over the 179.2-d interval spanned by the data. In the case of 
/s, the amplitude should increase from 3.9 to 9.6 mmag over 
this interval. A preliminary analysis of the y data showed 
that these predictions are in error, grossly so in the case of 
the first two modes. In order to examine this issue in de¬ 
tail, we divided the y, v, and u data into adjacent segments; 
except for the first segments in all Hlters and the last two 
segments in v and u, the segments partly overlapped. The 
hrst y segment spanned 40 d and each of the remaining y, 
V, and u segments, about 29 d. In each segment we then 
derived the amplitudes, Ai {i = 1, 2, 3), by the method of 
least squares using the following observational equations: 

N 

Arrij = Ao ^ Ai cos(27r/itj -I- (ki), (1) 

i = l 

where A = 3, Arrij are the y, v, or u differential magnitudes, 
fi {i = 1, 2, 3) are assumed to be equal to 5.9112, 5.8550 
and 5.5032 d~^, respectively, and tj are HJD ‘minus’ an ar¬ 
bitrary initial epoch. The assumed values of the frequencies 
fi {i = 1, 2, 3) are approximately equal to the weighted 
means of the values given in Tables 4, 5 and 8, respectively, 
with the weights equal to the corresponding amplitudes (see 
Section 7). Note that the hrst y segment covers about 2 cy¬ 
cles of the longest beat-period in the variation of 16 Lac, 
the one arising from the interference of /i and / 2 , and the 
remaining segments, about 1.5 cycle. An examination of the 
residuals from the least-squares solutions revealed six de¬ 
viant observations in the second y segment, all from the 
Mayaki Observatory. After rejecting these observations, the 
residuals in all but one segment fell within an interval of 
±0.015 mag; the exception was the hrst u segment in which 
there were hve residuals outside the ±0.015 mag interval, 
but all smaller in their absolute value than 0.018 mag. 
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Table 1. Log of the photometric measurements of 16 Lac. Observatories are listed in the order of their geographic longitude. 


Observatory 

Longitude 

Latitude 

Telescope 

Amount of data 
Nights h 

Filter (s) 

Observer(s) 

Sierra Nevada 

-3°23' 

37° 04' 

0.9-m 

18 

75.9 

uvby 

ER,PJA,RG 

Mercator 

-17° 53' 

28° 46' 

1.2-m 

33 

104.2 

Geneva 

KU,RD,JDD,TV, 

Fairborn 

-110° 42' 

31° 23' 

0.75-m APT 

57 

205.5 

uvy 

JDR,BA,POB 

Lowell 

-111°40' 

35° 12' 

0.5-m 

19 

96.8 

uvy 

MJ 

San Pedro Martir 

-115°28' 

31°03' 

1.5-m 

10 

36.4 

uvby 

EP,JPS,LP 

Sobaeksan Optical 

+128° 27' 

36° 56' 

0.6-m 

5 

31.9 

V 

SLK,JAL,SHK 

Astronomy 

Mt. Dushak-Erekdag 

+57° 53' 

37° 55' 

0.8-m 

12 

63.4 

V 

TND,NID 

Mayaki 

+30° 17' 

46° 24' 

0.5-m 

6 

13.0 

V 

AIM 

Piszkesteto 

+19° 54' 

47° 55' 

0.5-m 

13 

56.1 

V 

MP,DZ,DL,VA 

Bialkow 

+16°40' 

51° 29' 

0.6-m 

12 

67.2 

y 

ZK,GM,JMZ,AP,MS 

Total 




185 

749.3 
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Figure 1. A sample of our u, v, and y light-curves of 16 Lac (plus signs, diamonds, and circles, respectively). The lines are the synthetic 
light-curves computed as explained in Section 3. Note the deviation of the measurements from the fit due to three consecutive eclipses 
on the nights JD 245 2917, JD 245 2929, and JD 245 2941. The amount of data displayed is about half the total. 


Results of the above exercise are presented in Fig. [2] for 
the y and v amplitudes (points and circles, respectively), and 
in Fig. 13] for u. The error bars shown in the figures extend 
one standard deviation on each side of the plotted symbols. 
However, the standard deviations are not the formal stan¬ 
dard deviations of the least-squares solutions of equations 
(1) bu t twice the fo r mal s tanda rd deviations. In t h is, we 
follow iHandler et al.l ll200d) and jjerzvkiewicz et ^ (l2005l ) 


who—while dealing with time-series observations similar to 
the present ones—showed that the formal standard devia¬ 
tions were underestimated by a factor of about two. 

The solid line in the bottom panel of Fig. [2] shows the 
predicted blue-light amplitude of the /s mode (see the first 
paragraph of this section). The predicted blue-light ampli¬ 
tude of the fi mode, 12.4 mmag, is a factor of about 2 greater 
than the observed y and v amplitudes seen in the top panel 
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60 80 100 120 140 160 180 200 


HJD - 2452800 

Figure 2. The y (points) and v (circles) amplitudes of the /i, 
/ 2 , and /a modes of 16 Lac (from top to bottom) plotted as a 
function of HJD. The line in the bottom panel shows predicted 
blue-light amplitude of the /a mode (see the text for details). 



60 80 100 120 140 160 180 200 


HJD - 2452800 

Figure 3. The u amplitudes of the /i, / 2 , and /a modes of 16 
Lac (from top to bottom) plotted as a function of HJD. 


of Fig. 0 In the case of the /2 mode, the prediction fails on 
two accounts: (1) the observed amplitudes vary from uearly 
zero to about 0.8 of the predicted value of 15.1 mmag on 
a time scale about 50 times shorter than that determined 
by JP99, and (2) the median values of the observed y and 
V amplitudes are a factor of about 2 smaller than predicted 
blue-light amplitude. 


Standard frequency analysis and successive pre¬ 
whitening with sinusoids applied to the 2003-4 y magni¬ 
tudes yielded the following frequencies (in the order they 
were identified from the power spectra): 5.503, 5.591, 5.856 
and 5.852 d“^. The first three numbers are very nearly equal 
to /a, /i and / 2 , while the fourth is close to / 2 . The ampli¬ 
tudes amounted to 8.8, 5.9, 4.3 and 1.8 mmag, respectively. 
The first two amplitudes are close to the mean y amplitudes 
As and Ai, while the sum of the third and the fourth is 
close to the mean A 2 (see Fig. n. In addition, the fourth 
frequency differs from the third by less than the frequency 
resolution of the data. Clearly, the fourth frequency is an 
artefact. This shows that because of the variable amplitudes 
of the /2 and /a modes, the usual procedure of pre-whitening 
with sinusoids is inadequate in the present case. Therefore, 
instead of the usual procedure, we tried the following two 
methods: (1) pre-whitening separately in each segment, (2) 
pre-whitening with Ai assumed constant, A 2 assumed to 
vary quadratically with time, and A 3 , linearly. Both meth¬ 
ods were applied to the y data; in the case of v and u, we 
limited ourselves to method 1. In method 1, we computed 
residuals from the least-squares fit of equation (1) in each 
segment. Then, we took straight means of the residuals for 
a given tj in the overlapping parts of adjacent segments, so 
that for each tj there was one residual. In method 2, we 
used equation (1) with As = cxmst, A 2 = Bi + Citj + D\Aj, 
and As — B 2 C 2 tj. Since the first segment’s A 2 would 
not fit the quadratic relation derived from the remaining 
segments (see the middle panel of Fig. El), we applied the 
quadratic equation separately to this segment and then to 
the remaining data. In this case there was no need to aver¬ 
age residuals because the first and second segment do not 
overlap. Using the residuals as data, we then computed the 
amplitude spectra and the signal to noise ratio {S/N) as a 
function of frequency, where S is the amplitude for a given 
frequency, and N is the mean amplitude in 1 d“^ frequency 
intervals; in the first frequency interval, we omitted the am¬ 
plitudes for / l/T, where T is the total span of the data. 
For the y residuals, the signal-to-noise ratios are plotted in 
Fig. [4] as a function of frequency. We shall refer to the plots 
of this sort as the signal-to-noise spectra or S/N spectra. 
In Fig. 21 the S/N spectra are shown for method 1 and 2 
in the upper and lower panel, respectively. In both cases, 
the highest S/N peak occurs at 11.359 d“^, a frequency 
very nearly equal to the combination frequency f 2 + fs- The 
corresponding amplitudes amount to 0.56 and 0.59 mmag, 
and the S/N values are equal to 4.3 and 4.5, respectively. 
Thus, both peaks are significan t according to the popular 
criterion of iBreeer et al.l J 19931 ). In the S/N spectrum of 
method 1 residuals pre-whitened with /2 + fs, the highest 
peak, having the amplitude of 0.64 mmag and S/N = 4.2, 
occurred at 6.299 d“^. In the analogous spectrum of method 
2 residuals, there were no peaks with S/N > 4; the peak 
at 6.299 d“^ had S/N = 3.4. The frequency of 6.299 d“^ 
we shall refer to as /r. In the S/N spectrum of method 1 
residuals pre-whitened with the five frequencies, there was 
a peak at 0.085 d~^, rather close to /orb, but there was no 
peak at 2/orb, although one was present at this frequency 
in the power spec tra of the 1965 V data (see figure 1 of 
I Jerzvkiewic j fl99^ . The amplitude and S/N at /orb were 
equal to 0.55 mmag and 2.1, respectively. At 2/orb, the cor- 
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Figure 4. The signal-to-noise-ratio spectra, computed from the 
y data pre-whitened with the frequencies fi = 5.9112 d~^, /2 = 
5.8550 d~^, and /a = 5.5032 d~^ by means of method 1 (upper 
panel) and method 2 (lower panel). For details, see text. 


responding numbers were 0.34 mmag and 1.3, respectively; 
the 1965 V amplitude was equal to 2.1 mmag. 

In the case of the w-filter residuals, computed using 
method 1, the S/N spectrum showed the highest peak at 
6.301 d“^, very nearly equal to /r. In the second S/N spec¬ 
trum, obtained from the residuals computed with this fre¬ 
quency included in pre-whitening, the highest peak occurred 
at 11.359 d“^, the same combination frequency as that found 
in the y residuals. The S/N amounted to 3.8 and 3.9 in the 
first and the second spectrum , so that in the second spec¬ 
trum the iBreeer et al.l lll993l ) criterion of S/N > 3.5 for a 
combination frequency was satisfied. The amplitudes were 
now greater than the y-filter ones, viz. 0.61 and 0.64 mmag, 
respectively. At low frequencies, there were peaks close to 
/orb and 2/orb. The S/N and the amplitude at /orb amounted 
to 2.7 and 0.66 mmag, while at 2/orb, to 2.4 and 0.58 mmag. 
The phase of maximum light of the /orb term was 0.49 ± 0.04 
orbital phase, suggesting a reflection effect. In the case of the 
2/orb term, the phase of maximum light was 0.51 ± 0.04 or¬ 
bital phase, excluding an ellipsoidal variation as the cause. 
The 1965 B amplitude at 2/orb was equal to 1.6 mmag, and 
the phase of maximum light was 0.68 ± 0.02 orbital phase 
(see I Jerzvkiewica 1 19931 . table 7). In the case of the u-filter 
residuals, the highest peak in the first S/N spectrum was 
at 12.359 d“^, the -1-1 d“^ alias of the combination fre¬ 
quency /2 + f 3 - In the second S/N spectrum, the highest 
peak occurred at 1.092 d“^, and the second highest peak, 
at 6.301 d“^. At 11.359 d~^ in the first spectrum and at 
6.301 d“^ in the second spectrum, S/N were equal to 3.9 
and 3.1, respectively, and the /2 -I- /a and amplitudes 
were equal to 0.82 and 0.80 mmag. At /orb, the S/N and 
the amplitude amounted to 2.0 and 0.75 mmag, while at 
2/orb, to 1.0 and 0.38 mmag. 

The S/N spectrum of the y residuals, computed by 
means of method 1 but with all five significant terms (i.e. 
/i, / 2 , /a, /2 + /a, and /4) included showed no peaks higher 
than 3.8. We decided to terminate the frequency analysis 



-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 

Orbital phase 


Figure 5. The y, v and u residuals from the five-term pulsational 
solutions plotted as a function of orbital phase. The y deviant 
points (see Fig. n are not shown. 

at this stage. The y, v, and u fits computed with the five 
terms taken into account were used to plot the synthetic 
light-curves in Fig. [T] 

4 THE ECLIPSE 

4.1 The EBOP solutions 

The residuals, computed by means of method 1 with the 
five terms taken into account (see the last paragraph of the 
preceding section) for all y, v, and u observations, including 
those obtained during eclipses, are plotted in Fig. [5] as a 
function of orbital pha se. The ephemeris used in comp uting 
the phases was that of IPigulski &: Jerzvkiewica lll988l) . i.e. 

Min. light = HJD 243 9054.568 -F 12.09684E. (2) 

In Fig. [6] the y residuals are shown in a limited range of 
orbital phase around the primary eclipse (lower panel) and 
those around the phase of the secondary mid-eclipse, pre¬ 
dicted by the spectroscopic elements from solution IV of 
LOl (upper panel). No secondary eclipse can be detected: 
the mean y residual in the ±0.008 phase interval around 
the predicted mid-eclipse epoch amounts to 0.2 ± 0.4 mmag. 
Clearly, the secondary component is much fainter than the 
primary. 

In an attemp t to derive the parameters of 16 Lac we 
used lEtzel’sl (Il98ll ) computer program EBO P. The program, 
based on the Nelson-Davi s-Etzel model llNelson fc DavisI 
I 1972 I : IPopper fc Etzellll981^ . is well suited for dealing with 
detached systems such as the present one. The spectro¬ 
scopic parameters oj and e, needed to run the program, 
were taken from solution IV of LOl. The components were 
assumed to be spherical because the secondary compo¬ 
nent’s mass is much smaller than that of the primary and 
the system is well detached: for any reasonable assump- 
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tion about the primary component’s mass, the mass ra¬ 
tio would be equal to about 0.13 and the semimajor axis 
of the very nearly circular relative orbit, to about 50 Rq 
or 8 primary component’s radii. However, we included re¬ 
flected light from the secondary because a trace of a reflec¬ 
tion effect can be seen in Fig. [S] especially in y and v, in 
agreement with the results of the frequency analysis (see 
the fourth and the penultimate paragraphs of Section 3). 
The primary’s limb-darkening coefficients were interpolated 
from table 2 of Walter Van Hamm^il for TeS = 22 500 K, 
log q = 3.85 (see Section 5.1) and [M/H] = 0 llThoul et ^ 
l200d : iNiemczura fc Daszvhska-Daszkiewicj[^05ll . Unfortu- 
nately, since nothing is known about the secondary com¬ 
ponent except that it is much fainter than the primary, the 
central surface brightness of the secondary, Js, which EBOP 
uses as a fundamental parameter, must be derived indirectly. 
In units of the central surface brightness of the primary we 
have 


T _ 7-2 (1 Mp/3) la 

(1 - u,/3) /p ’ 


(3) 


where Ip and Is are the normalized lights of the primary and 
the secondary, respectively. Up and Ua are the limb darkening 
coefficients, and k = Vs/rp is the ratio of the radii. For y we 
have 


log= 2 log fc-I-4 log/ ct-I-ABC/2.5, (4) 

tp 

where fcx = Teff.s/Ieff.p is the ratio of the effective temper¬ 
atures of the components and ABC = BCs — BCp is the 
difference of the bolometric corrections. Introducing la/lp 
from this equation into equation (3) we get 

1 Us j 3 

The bolometric cor rection of the primary component was 
taken from table 3 of iFlowe 3 l|l99(^l for Teff.p = 22 500 K. As¬ 
suming Teff,s, we read the secondary’s bolometric correction 
from the same table as above. Assuming further loggr = 4.0 
and [M/H] = 0, we read the secondary’s limb-darkening coef¬ 
ficients from table 2 of Walter Van Hamme. Setting the inte¬ 
gration ring size to l°and the remaining parameters to their 
default EBOP values, we run the program for several values 
of Teft,s with Tp, the relative radius of the primary compo¬ 
nent, i, the inclination of the orbit, fc, the ratio of the radii, 
and Ss, the reflected light from the secondary as unknowns. 
As data, we used the y residuals (see Fig. [5] bottom panel). 
Unfortunately, we failed to find a solution which would con¬ 
verge. Convergent solutions were obtained if one of the first 
three unknowns, rp, i or fc, was fixed. After a number of tri¬ 
als we decided to fix fc, leaving rp, i, and Ss as the unknowns. 
For a given fc, identical triples of rp, i, and Ss were obtained 
for different T^fi^s, i.e. different Js- For example, for fc = 0.23, 
rp was equal to 0.1278 ± 0.0015 and i was 82?93±0?10, 
the same for Tes^s = 5000 and 6050 K; Sa was equal to 
0.00023 ±0.00005 for 5000 K, and to 0.00025 ± 0.00005 for 
6050 K. The synthetic light-curves computed from these so¬ 
lutions were very nearly identical everywhere but around 
the secondary eclipse: the computed depth of the secondary 


^ http://www2.fiu .edu/'^vanhamme/limdark.htm, see also 
IVan Hammj (Il993h 
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Figure 6. The y residuals from the five-term pulsational solution 
in a limited interval of orbital phase around the primary and 
secondary eclipse. The vertical line in the upper panel indicates 
the predicted epoch of the secondary mid-eclipse. The other lines 
are fragments of the synthetic light-curve, computed as described 
in the text. The open circles are deviant points not included in 
the analysis. 


eclipse was equal to 0.1 mmag for Teff.s = 5000 K, while 
it was 0.7 mmag for Tefi^s = 6050 K. The mean residual in 
the ±0.008 phase interval around the predicted mid-eclipse 
epoch amounted to 0.3 ± 0.4 and —0.2 ± 0.4 mmag for 5000 
and 6050 K, respectively. The mean residual was equal to 
0.0 ±0.4 mmag for 5750 K, and the computed depth of the 
secondary eclipse was then 0.5 mmag. For this computed 
depth of the secondary eclipse, the mean residual was equal 
to 0.0 ±0.4 mmag regardless of fc. 


4.2 Evolutionary state of the secondary 
component 

For a range of fc, synthetic light-curves computed from the 
solutions with an assumed depth of the secondary eclipse 
are indistinguishable from one another. Each solution yields 
an abscissa for plotting the secondary component in the HR 
diagram. For an assumed mass of the primary component, 
one can also have the secondary’s radius in absolute units 
from rp, i and the spectroscopic elements Ki and e, and 
therefore, the secondary’s ordinate in the HR diagram. For 
three values of the assumed depth of the secondary eclipse, 
0.1, 0.5 and 0.7 mmag, and the primary component’s mass 
Mp = 10 Mq, positions of the secondary component in the 
HR diagram are shown in Fig. [7] for a range of fc. 

For the range of fc shown in Fig. 0 the assump¬ 
tion of Mp = 10 Mq implies 1.30 ^ Ms < 1.31 Mq. 
Therefore, without making noticeable errors, we can com¬ 
pare the HR-diagram positions of the secondary compo¬ 
nent with M = 1.30 Mq evolutionary tracks. The evo¬ 
lutionary tracks plotted in the figure are the 1.30 Mq, 
Y = 0.26 5, Z — 0.0175 Pisa pre-main-sequence (pre- 
MS) tracks llTognelli. Prada Moroni fc Degrinnocentiir201ll . 
http://astro.df.unipi.it/stellar-models/). The tracks were 
computed using the mixing-length theory of convection with 
three values of the mixing length I = aHp, where a is the 
mixing-length parameter and Hp is the pressure scale height. 
The mid-value of Z (a = 1.68) was calibrated by means 
of the Pisa standard solar model (for further details, see 
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Figure 7. The HR-diagram positions of the secondary component 
of 16 Lac for three series of eclipse solutions, obtained with the 
ratio of the radii k = Vs/r-p as the parameter and assuming the 
depth of the secondary eclipse to be equal to 0.1 mmag (open 
circles), 0.5 mmag (dots) and 0.7 mmag (open squares). The open 
circles correspond to k from 0.21 to 0.31 with a step of 0.02, the 
dots, to k from 0.21 to 0.27 with a step of 0.02, and the open 
squares, to k equal to 0.21, 0.23, and 0.25; three open circles and 
one dot are labeled with their values of k. Also shown are Pisa 
1.3 Mq pre-MS evolutionary tracks for three values of the mixing- 
length parameter a (the solid line, the green dash-dotted line and 
the red dashed line). The thickened segment of the ol = 1.68 track 
indicates the age of 16.3 ±1.3 Myr, i.e. the evolutionary age of 
the primary component (see Section 5.2). 


iTognelli et al.|[201lh . The thickened segment of the a = 1.68 
track indicates the age within 1 ct of 16.3 Myr, the evolution¬ 
ary age of the primary component (see Section 5.2). Since 
the duration of the pre-MS phase of the primary compo¬ 
nent’s evolution is of the order of 0.1 Myr, we conclude 
that the secondary component is in the pre-MS contrac- 
tio n phase, the same conc l usion as that reached long ago 
bv IPigulski &: Jerzvkiewi^ ill9881 ). Now, the points of inter¬ 
section of the three series of the eclipse solutions in Fig. [7] 
with the pre-MS tracks constrain the ratio of the radii to 
0.23 < fc < 0.27. If the mass of the primary component were 
assumed to be equal to Mp = 8.8 Mq, corresponding to 
Ms « 1.20 Mq, the constraints would be 0.21 ^ k < 0.25. 
If Mp = 11.2 Mq (Ms « 1.40 Mq), 0.23 <k< 0.29. Thus, 
for 8.8 ^ Mp ^ 11.2 Mq we get 0.21 < fc < 0.29. Over this 
range of k, the relative radius of the primary component, 
rp, is a monotonically increasing function of k, while the 
inclination of the orbit, i, is monotonically decreasing with 
k, and both are virtually independent of the assumed depth 
of the secondary eclipse. Thus, from the last inequality we 
have 0.125 < rp < 0.132 and 83?4 > i > 82?0. More im¬ 
portantly, we can also obtain the lower and upper bound of 
the logarithmic surface gravity of the primary component: 
3-78 < loggp < 3.87. The formal standard deviations of 
rp, i and log(?p, equal to 0.0013, 0?10 and 0.011 dex, respec¬ 
tively, are—not surprisingly—much smaller than the allowed 
ranges of rp, i and loggp. 

The eclipse solutions which predict the depth of the 
secondary eclipse to be equal to 0.5 mmag yield synthetic 
light-curves which best ht the data around the secondary 
eclipse (see the end of the last paragraph of Section 4.1) 
while they ht the data elsewhere as well as do the other so¬ 
lutions. The fact that the corresponding line in Fig. [7] (the 
solid line with dots) crosses the a = 1.68 evolutionary track 


at an evolutionary age within the range of the evolutionary 
age of the primary component (see Section 5.2) is encourag¬ 
ing. It would be worthwhile to carry out space photometry 
of 16 Lac in order to hnd out whether the depth of the 
secondary eclipse is indeed close to 0.5 mmag. In any case, 
space photometry will be necessary to better constrain the 
range of k, and therefore the ranges of rp, i and loggfp. 


5 FUNDAMENTAL PARAMETERS 

5.1 The effective temperature and the surface 
gravity 


The effective temperature and surface gravity of 16 Lac 
can be obtained from the Stromgren indices using several 
photometric ca librations available i n the literature. The 
Cl index from iHauck fc MermilliodI (Il998l) . cor rected for 
the i nterstellar reddening in the standard way (ICrawfordI 
Il978ll . yielded the following values of Tbb (with the calibra- 
tio n referenced in the parenth eses after each value): 22 435 
K (iDavis fc Shobbroolj|l977^■ 22 580 K (UVBYBETAH). 
22 430 K dSterken fc Jerzvkiew^ 1993ll. ari d 22 635 K 
llBalonalll994l ). In the case of the lBalonal (119941 ) calibration, 
the P index was also needed in addition to cq. Taking a 
straight mean of these values we get = 22 520 K, with a 
formal standard error equal to 50 K. The latter number is so 
small because the four photometric calibrations are not in¬ 
dependent; th ey all rely h e avily on the OAO-2 absolute flux 
calibration of ICode et al.l (Il976l) . Realistic standard devia¬ 
tions of the effective temperatures of early-type stars, esti¬ 
mated from the uncert ainty of the absolute flux calibration, 
amou nt to about 3 % llNaniwotzki et al]ll993l : IJerzvkiewicj 
1 19941) or 680 K for the Teft in question. Thus, Teft of 16 Lac, 
obtained from the Stromgren indices, is equal to 22 520 ± 680 
K. The most recent spec troscopic determin ations of Tgfj 
incl ude 22 900 ±1000 K llThoul et al.l l2003l) . 21 500 ±750 
K-jPrugirrej^-Va ^l^^ Koleva 2011 ) and 23 000 ± 200 K 
llNieva fc Przvbillall2012l) . A straight mean of these numbers 
is equal to 22 470 ±480 K, in surprisingly good agreement 
with the photometric value. We shall adopt 22 500 ± 600 K 
as the Teft of 16 Lac. 

The logarithmic surface gravity of 16 Lac derived from 
P and Co turned out to be equal to 3.93 (UVBYBETA) 
and 3.90 llBalonalll994h . The good agreem ent between these 
values may be misleading: according to iNapiwotzki et aP 
lll993l) . the uncertainty of photometric surface gravities of 
hot stars is equal to 0.25 dex. We conclude that the photo¬ 
metric logg of 16 Lac is equal to 3.90 ±0.25. The spectro¬ 
scopi c values of log fir are equal to 3 .80 ± 0.20 llThoul et al.l 
l2003l). 3.75 ±0.17 i^gniel et all [ 2 ^ and 3.95 ±0.05 
I Nieva fc Przvbillal I 2 OI 2 I ) . A straight mean of these num¬ 
bers is equal to 3.83, with a standard error equal to 0.06 
dex. We shall adopt 3.85 ± 0.15 as the log g of 16 Lac, where 
the adopted standard deviation, equal to the median of the 


^ A FORTRAN prog ram _ based on the grid published 

by iMoon fc Dworetsk-^ lll985h . Written in 1985 by T.T. 
Moon of the University London and modified in 1992 
and 1997 by R. Napiwotzki of U niversitaet Kiel (see 
INapiwotzki. Schonberner fc Wenskelll993l) . 


© 0000 RAS, MNRAS 000, 000-000 



































































16 (EN) Lac and 2 And 9 



Figure 8. The primary component of 16 Lac plotted in the ef¬ 
fective temperature — surface gravity plane using Tefj and logg 
from Section 5.1 (dot with error bars). The dashed horizontal 
lines indicate the lower and upper bound of loggp, obtained in 
Section 4.2 from the eclipse solutions and pre-MS evolutionary 
tracks. The 9, 10 and 11 Mq evolutionary tracks are explained 
in the text. 

standard deviations of the individual log g values, is a com¬ 
promise between the standard deviation of the photometric 
log Q a nd that of the spectroscopic log g of iNieva &: Przvbillal 
(l2012t) . Note that because of negligible brightness of the sec¬ 
ondary component, the Tefi and logg we adopted pertain to 
the primary component. 


5.2 The effective temperature — surface gravity 
diagram 

In Fig. [S] the primary component of 16 Lac is plotted 
in the log Teg - logg plane using the effective tempera¬ 
ture and surface gravity from Section 5.1 (dot with error 
bars). The dashed horizontal lines indicate the lower and 
upper bound of loggfp, obtained in Section 4.2 from the 
eclipse solutions. Also shown are evolutionary tracks com¬ 
puted by mean s of the Warsaw - New Jersey evolutionary 
code (see, e.g.. IPamvatnvkh et al.l[l998l ). assuming the ini¬ 
tial abundance of hydrogen X = 0.7 a nd the metallicity Z 
= 0.0 15, the OPAL equation of state Jr^ rers &: NavfonovI 
I 2 OO 2 I ) and the OP opacities ([Seatonl 200^ for the latest 
heavy element mixture of lAsplund et all ll2009l ). For lower 
temperatures, the op acity data were supplemented with 
the Ferguson tables (iFerguson. Alexander fc AllardI l2005l : 
ISerenelli et al.ll2009l ). We assumed no convective-core over¬ 
shooting and Vrot = 20 km s“^ on the zero-age main 

sequence, a value consistent with t he ob served Dot sin i 
llGlebocki, Gnacihski fc Stawikowskil 120001 ) under the as¬ 
sumption that the rotation and orbital axes are aligned. 
The effect of rotation on log(? was taken into account by 
subtracting the centrifugal acceleration, amounting in the 
present case to about 0.001 dex. 

The evolutionary mass at the position of the dot in 
Fig. [8] Mev = 9.8± 1.3 Mq; the standard deviation of logg 
is responsible for the standard deviation of Mev. If the lower 
and upper bound of loggp were used, t he result would be 

9.3 5 Mev < 10.8 Mq. According to iNieva fc Przvbillal 
(l20ia) . Mev = 9.8±0.3 Mq. In this case, the small stan¬ 
dard deviation of Mev is the consequence of the small stan¬ 
dard deviations these authors assign to their spectroscopic 
Tefr and log;?. The evolutionary age of 16 Lac, obtained 


from Fig.[8l amount s to 16.3 ± 1.5 M yr, in surprisingly good 
agreement with the iBlaauw’d (1 19641 ) estimate of the age of 
Lac OB la. The most recent asteroseismic analysis of the 
star dThoul et al.ll200'^ ) has led to a mass of 9.62 ±0.11 Mq 
and an age of 15.7 Myr. Clearly, the accuracy of the astero¬ 
seismic values is much higher than that attainable by either 
the photometric or spectroscopic method. However, the sen¬ 
sitivity of the asteroseismic values to the details of modeling 
needs to be examined. We leave this for a future paper. 


6 THE HARMONIC DEGREE OF THE THREE 
HIGHEST-AMPLITUDE MODES 

6.1 From the uvy data 

Our data are the most extensive photometric observations 
of 16 Lac ever obtained. In addition, the three photomet¬ 
ric passbands we used include one on the short-wavelength 
side of the Balmer jump and two in the Paschen contin¬ 
uum. Thus, we can derive the amplitude ratios, AyjAu and 
A^/Au, and the phase differences, 4>j, — and "Fv — <I>„, 
which will be more accurate than any available before and 
sensitive to the harmonic degree of the pulsation modes. 
However, before computing the amplitudes and the phases 
we had to tackle the problem of the differences in the num¬ 
ber and the time distribution between the y, v and u data. 
Since the u data are fewer in number and less evenly dis¬ 
tributed in time than the y data, there were fewer u segments 
than the y segments (compare Fig. [3] with Fig. [^. Gonse- 
quently, the mean epochs of the u segments did not match 
those of the y segments. To a smaller degree, this was also 
the case with the v data. Because the amplitudes vary from 
one segment to another, the amplitude ratios computed us¬ 
ing amplitudes from unmatched segments would be biased; 
the same goes for the phases and the phase differences. We 
therefore divided the y and v data into new segments, in 
most cases different from those we formed in Section 3. In 
the new segments, the initial and final epochs and the time 
distribution of the data matched those of the u segments as 
closely as the data allowed. Then, the amplitudes and the 
phases in each segment were derived by fitting equation (1) 
with A = 5 to the data; the frequencies fi {i = 1,..,4) and 
/s = / 2±/3 were the same as in Section 3. Finally, the ampli¬ 
tudes of the three highest-amplitude modes from the match¬ 
ing segments were used to compute the amplitude ratios, 
and the phases, to compute the phase differences. In spite 
of the high quality of our photometry, standard deviations of 
the phase differences were rather large, rendering them use¬ 
less for mode identification. The amplitude ratios and their 
weighted means are listed in Table 2. There is no evidence 
in the table for a time-variability of the amplitude ratios, 
a result consistent with the fact that we are dealing with 
normal pulsation modes. In computing the weighted means, 
we assumed weights inversely proportional to the squares of 
the standard deviations of the components. Standard devi¬ 
ations of the weighted means were computed by adding the 
standard deviations of the components in quadrature and 
dividing the sum by the number of components. Note that 
for / 2 , the JD 245 2855.4 values are deviant. This is because 
at this epoch the /2 amplitudes were close to zero. In com¬ 
puting the weighted means, we omitted the JD 245 2855.4 
values. 
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Table 2. The uy and uv amplitude ratios for the three highest-amplitude modes. In the case of / 2 , the HJD 245 855.4 values were omitted 
in computing the weighted means. 


HJD-245 2800 

h 

Ay jA.U 

Av jAu 

/2 

Ay /Au 

Av fAu 

h 

Ay /Au 

Av j Au 

55.4 

0.565±0.078 

0.560±0.082 

0.324±0.870 

0.270±0.900 

0.6434i0.088 

0.6974i0.096 

107.1 

0.494±0.040 

0.540±0.040 

0.697±0.108 

0.843±0.114 

0.6494i0.040 

0.7234i0.040 

124.8 

0.489±0.028 

0.567±0.030 

0.753±0.066 

0.816±0.072 

0.6554i0.026 

0.7194i0.028 

135.1 

0.505±0.034 

0.588±0.040 

0.743±0.072 

0.8314:0.084 

0.6474i0.028 

0.7064i0.034 

160.0 

0.444±0.050 

0.494±0.056 

0.663±0.084 

0.6764:0.088 

0.6434i0.042 

0.6854i0.044 

198.4 

0.440±0.060 

0.517±0.070 

0.758±0.056 

0.8184:0.060 

0.6764i0.046 

0.7234i0.052 

Wt. Mean = 

0.488±0.021 

0.554±0.023 

0.734±0.035 

0.8004:0.038 

0.6524i0.020 

0.7124i0.022 




Figure 9. A comparison of the observed (symbols with error 
bars) and theoretical (lines) uvy amplitude ratios for the three 
highest-amplitude pulsation modes of 16 Lac, /i, /2 and f^. The 
observed amplitude ratios are the weighted means listed in the 
bottom line of Table 2. The theoretical amplitude ratios corre¬ 
spond to the dot in Fig. |8] The theoretical £ = 3 amplitude ratios 
are shown with the dashed line. 


Figure 10. A comparison of the observed (symbols with error 
bars) and theoretical (lines) UBG amplitude ratios for the three 
highest-amplitude pulsation modes of 16 Lac, /i, /2 and /a; the 
open square was shifted slightly along the abscissa to avoid over¬ 
lap. The observed amplitude ratios are the weighted means listed 
in the bottom line of Table 3. The theoretical amplitude ratios 
correspond to the dot in Fig. |8] 


A comparison of the observed amplitude ratios with the 
theoretical ones is presented in Fig. (9] The theoretical ampli¬ 
tudes were_corn£uted_aecording_tothe_zero-rotation formu¬ 
lae of iDaszviiska-Daszkiewic z et al.l (l2002l) usin g the nonadi- 
abatic pulsational code of iDziembowskil 1119771) and Ku rucz 
line-blanketed LTE model atmospheres ( Kuruc j l2004h for 
[M/H] = 0.0 and the microturbulent velocity ^ — 2 km/s. 
The remaining input parameters were the same as those 
used in computing the evolutionary tracks in Section 5.2. 
The calculations were carried out for log Tefj and log g used 
to plot the dot in Fig. |8] The loggfp obtained in Section 4.2, 
although more precise than logg, was not used because it 
may be less accurate on account of being model-dependent. 
As can be seen from Fig. |9l /i should be identihed with a 
radial mode, notwithstanding that the agreement between 
the observed and theoretical AyjAu is problematic. The re¬ 
maining two modes are nonradial with 1^3 because neither 
the ^ = 0, nor the £ = 4 line hts their amplitude ratios. In 
the case of /a, the observed and theoretical amplitude ratios 
agree to within 1 cr for ^ = 1, while in the case of / 2 , the 
observed A„/A„ falls about la below the 1 = 2 and 3 lines 
while the observed AyjA^ lies half-way between the i = 2 
and 3 lines. 


6.2 From the Geneva UBG data 

The data obtained with the Geneva filters (see Section 
2) cover three intervals: JD 245 2861.6 to JD 245 2872.7, 
JD 245 2921.5 to JD 245 2949.5, and JD 245 2971.3 to 
JD 245 2991.5. From the data in the latter two intervals we 
derived the amplitude ratios, Aq/Au and AsjAu in the 
same way as from the uvy data in Section 6.1. The ampli¬ 
tude ratios are listed in Table 3. The weighted mean am¬ 
plitude ratios are plotted as a function of the passbands’ 
central frequency in Fig. 1101 Also plotted are the theoreti¬ 
cal amplitude ratios. The harmonic-degree identification of 
£ = Q for /i and i = 1 for /a inferred in Section 6.1 from the 
uvy data is conhrmed. In the case of / 2 , £ = 2 hts now bet¬ 
ter than ^ = 3. In addition, 1 = 3 would be less satisfactory 
than 1 = 2 because of the effect of cancellation in integrat¬ 
ing over the stellar disc. However, the standard deviations 
of the /2 amplitude ratios are rather large. 

These harmonic-degree identihcations, i.e. £ = 0 for /i, 
£ = 2 OT, less satisfactorily, 1 = 3 for / 2 , and 1=1 for 
/a, agree with the earlier identihcations, based on the UBV 
amplitude ratios and the F-amplitude to the R V-amplitude 
ratio (see IPziembowski fc Jerzvkiewicj Il996h . They have 
points in common with the spectroscopic identihcations of 
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Table 3. The Geneva UG and UB amplitude ratios for the three highest-amplitude modes. 


HJD-245 2800 

h 

^gI-^u AbIAjj 

h 

AcjAu AbIAjj 

h 

AqIAu 

AbIAjj 

135.7 

0.477±0.086 

0.580±0.086 

0.727±0.155 

0.821±0.158 

0.647±0.060 

0.721±0.060 

181.9 

0.476±0.092 

0.563±0.090 

0.778±0.112 

0.833±0.108 

0.676±0.063 

0.717±0.060 

Wt. Mean = 

0.476±0.063 

0.572±0.062 

0.761±0.096 

0.829±0.096 

0.661±0.044 

0.719±0.043 


lAerts et al.! (l2003al l and lAerts et al.l ll2003bl '). An analysis 
of the line profiles of the He I A6678 A line led lAerts et al.l 
ll2003a^ to the conclusion that /i should be identified with 
a radial mode, / 2 , with an ^ = 2, m = 0 mode, and /a, with 
an ^ = 1, m = 0 mode. Subsequently. [Xerts et al] (l2003bl l 
modified the harmonic-deg ree identific a tion fo r fs to i < 3. 
For / 2 , the identification of lAerts et al.l (l2003bli is thus more 
specific than ours, while the reverse is true for /a. 


7 LONG-TERM VARIATION OF THE 

AMPLITUDE AND PHASE OF THE 

LARGE-AMPLITUDE TERMS 

7.1 The /i term 

All out-of-eclipse blue-filter observations of 16 Lac obtained 
throughout 1992 span an interval of over 40 years and con¬ 
sist of 6 334 data points (see JP99). By supplementing these 
data with our out-of-eclipse y observations we formed a data 
set of 9 384 points, spanning an interval of 53.4 years. We 
shall refer to this set as B&iy. Subtracting the contribution 
of the fi and fj terms (i,j = 1, 2, 3, i / j) from B&iy 
resulted in three sets which we shall refer to as R&j/ — 23, 
B&iy — 13 and B&iy — 12. We chose the archival blue-filter 
observations and the present y observations because they 
are much more numerous than observations in other filters. 
In the following frequency analysis of the combined data 
we shall neglect the difference between the blue and yel¬ 
low pulsation amplitudes. This will lead to some amplitude 
smearing in the results reported in this and the two fol¬ 
lowing sections. In 1965, when the amplitudes of the three 
high-amplitude terms were close to their maximum values, 
the difference between the B and V amplitudes amounted 
to 2.0±0.21, 0.40±0.2 1 and 0.70±0.20 m mag for /i, /2 
and /a, respectively fsee lJerzvkiewicjll993l . tables 7 and 8). 
Thus, the amplitude smearing will be the largest (albeit far 
from severe) in the case of fi and very nearly negligible in 
the remaining cases. 

The amplitude spectra of B&iy —23 are shown in Fig. 1111 
The abscissae of the highest peaks in the amplitude spec¬ 
tra are given in the caption to the figure and are listed in 
the second column of Table 4; the corresponding periods are 
given in column three. The amplitudes and phases with their 
formal standard deviations, obtained from a five-frequency 
least-squares fit of equation (1) to B&cy — 23 are listed in 
columns four and five. The standard deviation of the ht 
amounted to 5.4 mmag. The epochs of observations were 
reckoned from HJD 244 5784. 

The five-frequency fit accounts very well for the long¬ 
term variation of the amplitude and phase of the /i term. 
This can be seen from Fig. 1121 where the amplitude and the 
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Figure 11. The amplitude spectra of S&y—23, i.e. the 1950-1992 
archival blue-filter data and our ^-filter data combined, freed from 
the contribution of the /2 and /$ terms. Starting with the second 
panel from the top, the data were pre-whitened with (1) fi^i = 
5.911312 d-l, (2) /i,i and /i ,2 = 5.911275 d"!, (3) /i,i, /i ,2 and 
/ 1,3 = 5.911350 d-i, (4) /i,i, /i, 2 , /1,3 and /i ,4 = 5.911230 d'l, 
and (5) /i,i, /i, 2 , /l,3, /l ,4 and /i ,5 = 5.911395 d“l. 


Table 4. Fine structure of the fi term. 


3 

hj [d-'l 

Pl,j [d] 

Aij [mmag] 

4>lj [rad] 

1 

5.911312 

0.1691672 

12.45 ± 0.09 

0.181 ± 0.007 

2 

5.911275 

0.1691682 

9.06 ± 0.09 

4.080 ± 0.010 

3 

5.911350 

0.1691661 

3.27 ± 0.09 

3.326 ± 0.027 

4 

5.911230 

0.1691695 

2.26 ± 0.09 

1.836 ± 0.039 

5 

5.911395 

0.1691648 

1.35 ± 0.09 

5.367 ± 0.065 


phase of maximum light computed from the parameters of 
Table 4 (solid lines) are compared with the yearly mean 
amplitudes and the yearly mean phases of maximum light 
(upper and lower panel, respectively). The observed and 
computed phases of maximum light, ySmax, were obtained 
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Figure 12. Upper panel: The amplitude of the fi term, Ai, 
computed from the parameters of the five-frequency fit of Table 4 
(solid line), compared with the yearly mean blue-light 1950-1992 
Al (the dots, circles and the triangle) and the mean 2003 y and 
V Al (the filled and empty square, respectively). The dots, circles 
and the triangle are from figure 2 (top) of JP99 but with the error 
bars doubled. The 2003 error bars are not plotted because the}^ 
would be of about the same size as the symbols. Also shown is the 
amplitude, computed using the parameters of a fit with the first 
three frequencies of Table 4 (dashed line). Lower panel: The same 
for the phase of maximum light of the fi term. In both panels, 
the 2003 symbols are shifted along the abscissa by ±0.5 yr to 
avoid overlap. 


from the observed and computed epochs of maximum light, 
HJDmax, using the formula 

= 27v[E - (HJDmax - HJDq)/], (6) 

where E is the number of cycles which elapsed from an ar¬ 
bitrary initial epoch HJDo and / = /i,i from Table 4. The 
computed amplitude and phase of maximum light agree also 
with the nightly amplitudes and the nightly phases of max¬ 
imum light. This is illustrated in the upper half of Fig. 1131 
where the solid lines of Fig. [12] are plotted together with 
the 1965 B-filter amplitudes from JP96 and the 1965 B- 
filter phases of maximum light. The agreement between the 
computed amplitude and the y and v amplitudes derived 
in Section 3 is also satisfactory (see the upper panel of the 
lower half of Fig. o. The same goes for the computed and 
observed phases of maximum light (the lower panel of the 
lower half of the figure). Note that the lines in Figs. 1121 and 
were not fitted to the points shown in the figures but 
were computed independently from the parameters of the 
five-frequency fit. From Fig. [T^jit is also clear that the three 
first frequencies of Table 4 alone would be insufficient to ac¬ 
count for the variation of the phase, especially around 1900. 
This is an important conclusion because the first three pe¬ 
riods in Table 4, Pi,i, Pi,2 and Pi,3, are very nearly equal 
to the periods Pi, and of LOl, mentioned in the 
Introduction: the differences (in the sense ‘Table 4 minus 
LOT) amount to 0.00000013, 0.00000011 and 0.00000005 d, 
respectively. 

The first three frequencies of Table 4 form a very nearly 
equally-spaced triplet, /i,2, /i,i, /i,3, with a mean spacing 
equal to 0.0137 yr“^. The remaining frequencies, /i,4 and 
/i,s, flank the triplet at a distance of 0.0164 yr“^ from the 
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Figure 13. Upper half, upper panel: The amplitude of the /i 
term, Ai, computed from the parameters of the five-frequency 
fit of Table 4 (solid line), compared with the nightly 1965 B- 
filter Al from JP96 (dots). Lower half, upper panel: The same 
for the y and v Ai derived in Section 3 (the dots and circles, 
respectively). Upper half, lower panel: The phase of maximum 
light of the fi term, computed from the parameters of the five- 
frequency fit of Table 4 (solid line), compared with the 1965 B- 
filter phases of maximum light (dots). Lower half, lower panel: 
The same for the y and v phases of maximum light (the dots and 
circles, respectively). 

first and the last frequency of the triplet, respectively. The 
triplet’s spacing implies a time-scale of 73 yr for the long¬ 
term variation of the amplitude and phase of the /i term, 
while accounting for the non-sinusoidal shape of the varia¬ 
tion seen in Fig. 1121 requires all five components. Note that 
the reciprocal of the time-span of the data is equal to 1/53.4 
= 0.0187 yr“^, so that the spacings of the adjacent frequen¬ 
cies of the quintuplet amount to about 3/4 of the formal fre¬ 
quency resolution of the data. That they could be resolved 
nevertheless is due to their unequal amplitudes. 

7.2 The /2 term 

Frequency analysis of B&iy — 13 yielded four frequencies be¬ 
fore the noise prevented detecting further frequencies. How¬ 
ever, a four-frequency fit did not account very well for the 
variation of the amplitude and phase of the /2 term. After 
a number of trials, we found that the fit improved when 
the amplitudes were assumed to vary uniformly with time, 
i.e. when constant amplitudes in the observational equations 
were replaced by A 2 ,j + B 2 ,jt. The parameters of a least- 
squares fit with the amplitudes modified in this way are 
listed in Table 5. The standard deviation of the fit amounted 
to 5.6 mmag. The epochs of observations were reckoned from 
HJD 244 5784. A comparison of the amplitudes and phases 
of maximum light computed from the parameters of Table 
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Figure 14. Upper panel: The amplitude of the /2 term, A 2 , 
computed from the parameters of the four-frequency fit of Table 
5 (solid line), compared with the yearly mean blue-light 1950-1992 
A 2 (the dots, circles and the triangle) and the mean 2003 y and 
V A2 (the filled and empty square, respectively). The dots, circles 
and the triangle are from figure 3 (top) of JP99 but with the error 
bars doubled. The 2003 error bars are not plotted because they 
would be of about the same size as the symbols. Lower panel: 
The same for the phase of maximum light of the /2 term. In 
both panels, the 2003 symbols were shifted along the abscissa by 
±0.5 yr to avoid overlap. 

5 (solid lines) with the observed amplitudes and phases of 
maximum light is shown in Figs. [14] and 1151 The observed 
and computed phases of maximum light were obtained from 
the observed and computed maxima using equation (5) with 
/ = /2,i from Table 5. The agreement between the computed 
and observed amplitudes and phases of maximum light of 
the /2 term seen in Figs. I14|and ll5l is less satisfactory than 
was the case for /i (Figs. fT^ and [T 3 ]) . 

As can be seen from Table 5, /2,i — f 2,2 ~ 0.002607 
and /2,4 — /2.1 = 0.002648 d“^, so that /2,2, /2,i and /2,4 
form a very nearly equidistant frequency triplet, with a mean 
separation equal to 0.002628 d“^. Although this number is 
close to 1 yr“^, it is not an artefact because the aliases were 
removed in our procedure of pre-whitening. The beat period 
corresponding to the mean separation of the triplet is equal 
to 380.5 d, a value about 10 % greater than the beat-period 
between the LOl periods P 2 and . The difference between 
/ 2,2 and / 2,3 is equal to 0.000063 d“^ or 0.023 yr“^, implying 
a time scale of 43 yr. The latter number is close to the time 
scale of the variation of A 2 derived by JP99 from the 1950- 
1992 data. 

7.3 The /a term 

As we mentioned in the Introduction, the fs term was found 
by JP96 to be a doublet. Using all blue-filter observations 
of 16 Lac available at the time, JP99 determined the dou¬ 
blet frequencies to be /a,! = 5.5025779 ± 0.0000005 and 
/a,2 = 5.5040531 ±0.0000008 d“^. Slightly different frequen¬ 
cies, equal to 5.5025928 and 5.5040765 d“^, were subse¬ 
quently derived by LOl from RV data. A frequency analysis 
of BSzy — 12 showed this term to be a triplet. Using the 
frequencies read off the amplitude spectra as starting val¬ 
ues in a three-frequency nonlinear least-squares fit of equa¬ 
tion (1) to B&iy — 12 resulted in the frequencies, ampli¬ 
tudes and phases listed in Table 5. The standard deviation 
of the fit amounted to 5.4 mmag. The epochs of observations 
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Figure 15. Upper half, upper panel: The amplitude of the /2 
term, A 2 , computed from the parameters of the four-frequency 
fit of Table 5 (solid line), compared with the nightly 1965 B-filter 
A 2 from JP96 (dots). Lower half, upper panel: The same for the y 
and V A 2 derived in Section 3 (the dots and circles, respectively). 
Upper half, lower panel: The phase of maximum light of the /2 
term, computed from the parameters of the five-frequency fit of 
Table 5 (solid line), compared with the 1965 B-filter phases of 
maximum light (dots). Lower half, lower panel: The same for the 
y and v phases of maximum light (the dots and circles, respec¬ 
tively). The first-segment y and v phases of maximum light, equal 
to —1.6 ± 1.0 and 1.7 ± 1.7 rad, respectively, are not shown. 


were reckoned from HJD 244 5784. A comparison of the am¬ 
plitudes and phases of maximum light computed from the 
parameters listed in Table 5 with observations is shown in 
Figs. ll6l and ll7l The observed and computed phases of max¬ 
imum light were obtained from the observed and computed 
maxima using equation (5) with / = /a,i from Table 5. 

The new values of /a,i and /a, 2 are very nearly equal to 
those obtained by JP99. The new value of the beat-period is 
equal to 674.94 ±0.30 d, in good agreement with the beat- 
period between the LOl’s periods Pa and P^, mentioned in 
the Introduction. The third frequency, placed asymmetri¬ 
cally between /a,i and /a, 2, gives rise to two beat-periods, 
720.73 ±0.35 d and 29.09 ±0.43 yr. 
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Table 5. Fine structure of the /2 term. 


j 

/ 2 ,i [d- 1 ] 

P 2 J [d] 

A 2 J [mmag] 

B 2 J [mmag/d] 

^ 2 ,j [rad] 

1 

5.855574 

0.1707775 

6.828 ± 0.093 

-0.000173 ± 0.000021 

5.017 ± 0.013 

2 

5.852967 

0.1708535 

3.592 ± 0.096 

0.000047 ± 0.000016 

5.199 ± 0.028 

3 

5.852904 

0.1708554 

3.008 ± 0.090 

-0.000074 ± 0.000023 

5.483 ± 0.028 

4 

5.858222 

0.1707003 

2.164 ± 0.096 

0.000108 ± 0.000016 

3.836 ± 0.043 


Table 6 . Fine structure of the /a term. 


j 

/3,, [d- 1 ] 

Paj [d] 

Asj [mmag] 

[rad] 

1 

5.50257795 ± 0.00000033 

0.181733000 ± 0.000000011 

7.97 ± 0.10 

1.914 ± 0.013 

2 

5.50405957 ± 0.00000058 

0.181684080 ± 0.000000019 

4.24 ± 0.10 

2.484 ± 0.024 

3 

5.50396544 ± 0.00000127 

0.181687187 ± 0.000000042 

1.84 ± 0.09 

2.348 ± 0.049 



Figure 16. Upper panel: The amplitude of the /a term, A 3 , 
computed from the parameters of the three-frequency fit of Table 
5 (solid line), compared with the yearly mean blue-light 1950- 
1992 A 3 from table Al of JP96 and table 1 of JP99 (the dots and 
circles; the latter represent data of lower weight) and the mean 
2003 y and v A 3 (the filled and empty square, respectively). The 
2003 error bars are not plotted because they would be of about the 
same size as the symbols. Lower panel: The same for the phase of 
maximum light of the /a term. In both panels, the 2003 symbols 
were shifted along the abscissa by ±0.5 yr to avoid overlap. 


8 SUMMARY AND DISCUSSION 

Over the 179.2-d interval spanned by the present multisite 
uvy data, the amplitude of the fi term was constant but the 
amplitudes of the /2 and /a terms, A 2 and A3, increased by 
several mmags (see Figs. [2] and [3]). The latter fact made the 
usual procedure of pre-whitening inapplicable. In Section 3, 
using the values of the three frequencies determined from 
all available B, b and y data (Section 7), we pre-whitened 
the data piecewise, in segments so long that the three terms 
could be resolved but short enough to neglect the variation 
of A2 and A3. In the pre-whitened data we detected two 
low-amplitude terms having frequencies = 11.3582 and 
/s = 6.2990 d“^. The former frequency is equal to /2 ± /a, 
the latter is new. The amplitudes of the term amount to 
0.57 ± 0.09, 0.64 ± 0.10 and 0.92 ± 0.15 mmag for y, v and u, 
respectively, while those of /s, to 0.64 ± 0.09, 0.61 ± 0.10 and 
1.05 ±0.15 mmag, respectively. The frequencies 2/orb, 2/i, 
fi ± /s and fi ± /2 seen in 1965 (see the Introduction) were 
not found. Why /2 ± /s was present in 1965 and 2003-2004 
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Figure 17. Upper half, upper panel: The amplitude of the /a 
term, A 3 , computed from the parameters of the three-frequency 
fit of Table 5 (solid line), compared with the nightly 1965 B- 
filter Aa from JP96 (dots). Lower half, upper panel: The same 
for the y and v A 3 derived in Section 3 (the dots and circles, 
respectively). Upper half, lower panel: The phase of maximum 
light of the /a term, computed from the parameters of the three- 
frequency fit of Table 5 (solid line), compared with the 1965 B- 
filter phases of maximum light (dots). Lower half, lower panel: 
The same for the y and v phases of maximum light (the dots and 
circles, respectively). 
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while 2/i, /i + /a and /i + /2 were absent in 2003-2004 is 
easy to understand: in 2003-2004 A 2 was moderately smaller 
and As was slightly smaller than in 1965 (see Figs. 1141 and 
1161) while Ai decreased between 1965 and 2003-2004 to about 
one fourth of its 1965 value (see Fig. 11211 . The fact that 2/orb 
was missing in 2003-2004 suggests that it may be related to 
/i. Finally, /s was not detected in 1965 because its V and 
B amplitudes were then below ~0.40 mmag. 

With the fi {i = 1,.., 5) terms taken into account, we re¬ 
peated the piecewise pre-whitening. The resulting residuals 
were used in Section 4 to plot the eclipse light-curves (Figs. [5] 
and 01. The light curves show no ellipticity effect and no 
secondary eclipse can be detected. However, a marginal re¬ 
flection effect is present. From the y light-curves and the LOl 
spectroscopic elements uj and e, we computed the relative ra¬ 
dius of the primary, rp, the orbital inclination, i, and the re¬ 
flected light from the secondary^ 5's, for a range of fe = r^/vp 
by means of EBOP (lEtzellll98lh . This allowed an examina¬ 
tion of the position of the secondary component in the HR 
diagram in relation to pre-MS evolutionary tracks (Fig. 0, 
leading to the conclusions that (1) the secondary component 
is in the pre-MS phase of its evolution, and (2) the param¬ 
eters of the system can be constrained to 0.21 < k < 0.29, 
0.125 < rp < 0.132 and 83?4 > i > 82?0. 

For /i, /2 and /a, the uvy and the Geneva UBG ampli¬ 
tude ratios are derived from the multisite data and compared 
with the theoretical ones for the spherical-harmonic degree 
£ = 0,..,4 in Section 6. The theoretical amplitude ratios were 
computed using Teff and log g of 16 Lac obtained in Section 
5. In Section 6, the highest degree, £ = 4, is shown to be in¬ 
compatible with the observations. The first term, /i, could 
be identified with an £ = 0 mode, while the third, /a, with 
an £ = \ mode. In the case of /2, an unambiguous spherical- 
harmonic degree identification was not possible: it can be 
either an ^ = 2 or 3 mode, with the latter possibility less 
likely because of the effect of cancellation in integrating over 
the stellar disc. 

In Section 7, using the present y-HIter magnitudes and 
archival blue-filter magnitudes, we investigate the long-term 
variation of the amplitudes and phases of the three high- 
amplitude terms over the interval of 53.4 yr spanned by the 
data. In the case of /i, the magnitudes can be represented 
by means of a sum of five sinusoidal components with closely 
spaced frequencies (see Table 4). The first three frequencies 
form an equally-spaced triplet with a spacing of 0.0137 yr“^, 
implying a time-scale of 73 yr, in agreement with JP99 and 
LOl. The sum of the five components accounts very well for 
the non-sinusoidal shape of the variation seen in Fig. 1121 
Since the fi mode is radial, the regularities in the frequency 
spacings suggest an underlying amplitude and phase modu¬ 
lation of a single pulsation mode. 

In the case of /2 we could represent the 1950-2003 mag¬ 
nitudes by a sum of four sinusoidal components with uni¬ 
formly variable amplitudes (see Table 5). Of these, three 
components (/2,2, /2,i and /2,4 in the order of increasing fre¬ 
quency) form a triplet very nearly equidistant in frequency 
with a mean separation of 0.002628 d“^, corresponding to 
a beat period of 380.5 d, while the /2,3 component precedes 
the first member of the triplet by 0.023 yr~^, corresponding 
to a beat-period of 43 yr. The 380.5-d beat-period dominates 
the variation of the amplitude and phase (see Fig. [Til) : the 
43-yr beat-period is close to the time scale of the variation 


of A 2 derived by JP99 from the 1950-1992 data. The /2,2, 
/2,i, /2,4 triplet may be the result of (1) rotational splitting 
of a nonradial mode, (2) an accidental coincidence of two 
nonradial modes, one rotationally split, and (3) an acciden¬ 
tal coincidence of three nonradial modes. In all cases, only 
some members of the rotationally split multiplets would be 
excited to observable amplitudes. In (1) and (2), the veloc¬ 
ity of the star’s rotation can be computed from the Ledoux 
first-order formula using the radius from the eclipse solu¬ 
tions. The results are Kot < 1 km s“^ if (1), and Kot 2 
km s“^ if (2). Thes e numbers are in seve re conflict with the 
observed Frotsini (iGlebocki et ahlbood l. Thus, we are left 
with (3) and the problem why three accidentally coincident 
modes should be nearly equidistant in frequency. Finally, 
there is the possibility that our representation by means of 
a sum of four sinusoidal components may be merely a for¬ 
mal description of the complex long-term variation of the 
amplitude and phase of the /2 term. 

JP96 have noted that the reciprocal of the growth rates 
of several £ ^ 2 pulsation modes which are excited in models 
of 16 Lac and can be identified with the /i and /2 terms are 
of the same order of magnitude as the time scales of 73 and 
43 yr just mentioned. Since amplitude and phase modulation 
on a time-scale of the order of the reciprocal of the growth 
rate is predicted by the theory of non-linear interaction of 
pulsation modes, JP96 suggested that these time scales in 
the observed long-term behaviour of the /i and /2 modes 
result from (1) the 1:1 resonance between them, or (2) a 
resonant coupling to other modes. These points still await 
theoretical verihcation. 

Finally, the /a 1950-2003 magnitudes could be repre¬ 
sented by a sum of three sinusoidal components (see Table 
5). The hrst two components have frequencies /s.i and /a,2 
very nearly equal to those derived by JP96 and JP99 from 
the data available at the time. These two frequencies pro¬ 
duce a beat period of 720.73 ±0.35 d which dominates the 
amplitude and phase variation (see Fig. Hi). The frequencies 
/a,a and /a, 2 produce a beat period of 29.09 ± 0.43 yr. As in 
the case of /2, assuming that the /a,i, /a, 2 doublet is due 
to rotational splitting of a nonradial mode leads to Kot ^ 1 
km s“^, much smaller than observed. JP96 suggested that 
the doublet represents an accidental near-coincidence of two 
self-excited modes. In view of the unambiguous identihca- 
tion of the harmonic degree of the /a term (Section 6), these 
two modes must both have £ = 1. On the other hand, the 
members of the /2,2, /2,i, /2,4 triplet may have different £, 
a situ ation encoun tered earlier in 1 M on (iBalona fc Stobia 
I198OII and 12 Lac ([Handler et al.l[2006l L 
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APPENDIX A: THE VARIABLE COMPARISON 
STAR 2 ANDROMEDAE 

Al Frequency analysis 

As mentioned in the Introduction, the comparison star 
2 And turned out to be a small-amplitude variable. Fig. lAll 
shows the results of a frequency analysis of the differential y 
magnitudes 2 And — 10 Lac. The less accurate Sierra Nevada 
and Piszkesteto observations were omitted. 

The variability of 2 And can be described by a sin¬ 
gle frequency (Fig. lAll) . which is confirmed by the anal¬ 
ysis of the measurements in the other Alters. Nonlinear 
least-squares fitting to our j/-fllter data gives a value of 
7.195301 ±0.00012 d“^, whereby the error estimate was 
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Figure Al. Amplitude spectrum of 2 And. The top panel shows 
the spectral window of the data, followed by the periodogram of 
the data in the middle panel. Pre-whitening with /i leaves only 
noise in the residuals (bottom panel), aside from some residual 
extinction effects (peaks at 2.0, 3.0, 4.0 and 5.0 d~^). 


Table Al. Amplitudes and relative phases of the single frequency 
7.195301 d“^ for the var iable comparison star 2 And. For mal 
error estimates (following iMontgomerv 0’Donoghuelll999r) for 
the amplitudes and phases are given. 


Filter 

Amplitude 

[mmag] 

<j>x — (f>y 

[rad] 

u 

1.72±0.13 

-0.16 ±0.09 

V 

1.13±0.08 

-0.05 ±0.08 

y 

1.91 ±0.07 

0 by definition 

u 

2.07 ±0.29 

0.07±0.14 

Bl 

1.57±0.25 

-0.17±0.16 

B 

1.69 ±0.25 

-0.06 ±0.15 

B2 

1.65 ±0.26 

-0.22 ±0.16 

VI 

2.26 ±0.25 

-0.02 ±0.12 

V 

2.05 ±0.22 

0.01±0.12 

G 

2.42 ±0.26 

-0.04±0.11 


compu ted with the formulae of I Montgomery Sz Q’Donoghud 
(| 19991) . Fitting a sine curve of this frequency to the data re¬ 
sults in the amplitudes and relative phases, (px — (pv {x = u, 
V, U, Bl, B, B2, VI, V, G), listed in Table Al. 


A2 The cause of variability of 2 And 

Given the spectral type of A3 Vn and the period of the vari¬ 
ability, one would immediately suspect that the light vari¬ 
ations of 2 And are due to 5 Scuti-type pulsatio n. In fact, 
the star was so classihed by [Handler et all (l2006l l. However, 
if 2 And were a S Scuti variable, the amplitudes are gen¬ 
erally expected to increase from y to v and from G to Bl, 
and then level out or drop again towards u and U (e.g.. 


see lHevnderickx|ll994h . Interestingly, just the inverse is the 
case. In the following we investigate possible causes for this 
observation. 

2 And is not a single star. It is a close vis ual binary 
discovered over a century ago (iBurnhamI 1 18941 . The two 
components are physically associated. The latest determi¬ 
nation of the orbital period is 74 yr with an eccentricit y 
of 0.8 and a semimajor axis of 0.23” (H ica Romerdl2010l ~). 
Transf orming the Tycho-2 photo metry of the two compo¬ 
nents (iFabricius fc Makarov! 2000lj. to the standard Johnson 
system according to iBesseOI I 2000ll gives; Va = 5.24 mag, 
{B — V)a = 0.07 mag for 2 And A, and Vb = 7.51 mag, 
{B — V)b = 0.23 mag for 2 And B. Therefore Va-i-b = 5.113 
mag and [B — H)a+b = 0.086 mag, in reasonable agreement 
with the measu red V = 5.100 ma g and B — V = 0.094 mag 
for the system llMermilliodlll99ll j. 

_ The revised Hipparcos parallax of 2 And dvan LeeuwerJ 

l2007h implies a distance of 129 ± 9 pc. Concerning redden¬ 
ing, a comparison o f the results from the galactic reddening 
law of IChen et al.l (Il998l l for 2 And, and the reddening of 
two stars within 4° of 2 And in the sky, at a similar Hip¬ 
parcos distance (HR 8870 and HD 218394) lead us to adopt 
E{B-V) = 0.022 ±0.022 mag and A(H) = 0.07 ±0.07 mag. 

We then arrive at My = —0.39 ±0.16 mag for 2 
And A and 1.88 ±0.16 mag for 2 And B. The relations of 
I Flower! (Il996lj yield estimates of Teff = 8950 ± 250 K and 
Mboi = —0.45 ± 0.16 mag for 2 And A, and 7720 ± 250 K and 
1.85 ± 0.16 mag for 2 And B. The positions of the two com¬ 
ponents are shown in a theoretical HR diagram in Fig. IA2I 
The evolutionary tracks plotted in this hgure were com¬ 
puted with the Warsaw-New Jersey stellar evolution code, 
the OPAL equation of state and the OPAL opacity tables 
lllglesias fc Roger^ll996h , a hydrogen abundance of X = 0.7 , 
a metal abundance ot Z — 0.012 llAsplund et al.l l2004l L 
and no convective core overshooting. We assumed Vjot = 
250 km s“^ on the zero-age main sequence, so that the ob - 
served Kot sin i = 212 km s”*^ llRover. Zorec fc GomedlioO^ 
would be approximately matched. 

One finds that 2 And A is an M = 2.7 ±0.1 Mq star 
that probably has already left the main sequence (logg = 
3.40 ± 0.12, log age = 8.57 ± 0.04). Including core overshoot¬ 
ing would decrease M and the age but logg would remain 
virtually unchanged. 2 And A is located off the 5 Scuti 
instability strip. 2 And B is situated in the centre of the 
strip (M = 1.78 ±0.06 Mq, loggr = 3.90 ±0.16, log age = 


8.84 


+0.07 


With these parameters in hand, we shall test four pos¬ 
sible hypotheses to explain the cause of the variability of 2 
And: 5 Scuti pulsation or ellipsoidal variability, of either 2 
And A or 2 And B. Rotational modulation can be excluded 
immediately because either star would have to rotate faster 
than its breakup velocity. 

Regarding the hypothesis of <5 Scuti pulsation, we com¬ 
puted theoretical pulsation amplitudes of 5 Scuti models 
for the Stromgren and Geneva passbands in the parameter 
ranges depicted in Fig. lA2l for an ass umed radial veloc i ty am - 
plitude, using the methodology bv iBalona fc Evers! (Il999!) . 
We note that whilst computing realistic absolute amplitudes 
is still out of reach, amplitude ratios are easily obtainable. 

We computed theoretical pulsation amplitudes of 
S Scuti models for the Stromgren and Geneva passbands. 
From the flux ratio of 2 And A to 2 And B in the y band (es- 
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Figure A2. The positions of 2 And A and B in the theoreti¬ 
cal HR diagram. Some stellar model evolutionary tracks labeled 
with their masses in solar units (solid lines) are i ncluded. The bor¬ 
ders of the 5 Scuti instability strip (adopted from iRodrfguez et al.l 
I 2 OOCII . dotted lines) are included for comparison. The solid diago¬ 
nal line is the ZAMS. 


timated from Va and Vb) and the observed y amplitude from 
Table Al we then determined the undiluted photometric am¬ 
plitudes of 2 And A and 2 And B in the Stromgren y band, 
(Atrue = 0.0022 and 0.0215 mag, respectively), and used the 
model amplitudes to predict the intrinsic amplitudes in the 
other bands. As the next step, we determined the flux ratio 
of the two compone nts of 2 And in the different passbands. 
To this end, we used iKurucd (1 19941 ) model atmospheres, rep¬ 
resenting 2 And A with a Teg = 9000 K, logg = 3.5 model 
atmosphere, and 2 And B with a Teff = 7750 K, log 5 = 4.0 
model atmosphere. We then integrated the monochromatic 
fluxes from these model atmospheres over the ten photo¬ 
metric passbands used, and scaled the resulting fluxes to 
the observed ratio in y. 

The results of these computations can be summarized 
as follows: the behaviour of the pulsation amplitudes with 
wavelength is inconsistent with the observations for any kind 
of assumed pulsation of 2 And A. The reason is that the 
temperature of the star is so high, that the pulsation ampli¬ 
tudes always increase towards shorter wavelength, no mat¬ 
ter which spherical degree of the pulsation was tested (we 
checked for 0 ^ ^ 8). The additional flux of 2 And B does 

not change this behaviour due to the considerable luminosity 
difference. 

Assuming that 2 And B was a 5 Scuti pulsator, and 
considering the amplitude contamination due to the light 
of 2 And A, we find that the observed pulsation amplitude 
vs. wavelength dependence is only consistent with modes 
of even spherical degree of f ^ 6. We recall that in such 
a scenario the intrinsic pulsation amplitude of 2 And B 
in the Stromgren y band, Atme is 0.0215 mag. Geomet¬ 
ric cancellation of modes with I is 6 decreases the ob- 


served amplitudes to less than 1/5 0 of the intrinsic value 
lIPaszvhska-Daszkiewicz et al.ll20^ ). i.e. the intrinsic pulsa¬ 
tion amplitude would be enormous in this case. We therefore 
consider <5 Scuti pulsation of either of the two A-type stars 
in the 2 And system as highly unlikely. 

A remaining possibility is an ellipsoidal variation of ei¬ 
ther component. Consequently, we repeated the previous 
procedure under the assumption of no colour dependence 
of the amplitude. We found that the resulting wavelength- 
amplitude dependence explains the observations much bet¬ 
ter than the 5 Scuti-pulsation hypothesis. This implies an 
orbital period of twice the observed value. Using the effec¬ 
tive temperatures, the luminosities and the masses derived 
above we find from the Kepler’s third law that in the case of 
2 And A the sum of the radii of the components Aa and Ab, 
7?Aa + RAh, would be a factor of about 2 greater than the 
semimajor axis of the relative orbit, a, regardless of whether 
we assume (1) equal luminosity and mass components, or 
(2) the secondary component much fainter and less massive 
than the primary component. In the case of 2 And B, as¬ 
sumption (1) leads to i^Ba + R^h ~ 1.5o, while assuming 
(2) one gets i?Ba ~ 1.0a. Thus, an ellipsoidal variation of 2 
And B, arising from a tidal distortion of component Ba by a 
much less massive secondary component Bb in a tight orbit 
may be the cause of the variability of 2 And. 

In order to examine the last possibility let us use 
the lowest order approximation, certainly adequate in the 
present case, according to which the amplitude of the light 
variation due to aspect changes of the tidally distorted pri¬ 
mary can be expressed by means of the following formula: 

5m — 1.629AAg(i?Ba/a)® sin^ i, (Al) 

where Sm is expressed in magnitudes, A\ is the p hotomet- 
ric distortion parameter of iRussell fc Merrilll 1119521) . q is the 
mass ratio of the components, and i is the inclination of the 
orbit to the tangent plane of the sk y. We note tha t this for¬ 
mula is equivalent to equation (6) of iMorri i (1 19851 ). Neglect¬ 
ing the light variation of the secondary, we have 5m = Atrue. 
Then, taking into account the fact that in the present case 
A\{R^s./a)^ « 1, we get q « 0.013 and 0.052 for i = 90°and 
30°, respectively. The corresponding masses of the secondary 
component are equal to >^0.03 and ~0.11 Mq, i.e. they are 
in the range of the masses of brown dwarfs. Thus, our hy¬ 
pothesis that 2 And B is an ellipsoidal variable leads to the 
following model: the 2 And Bb component is a brown dwarf 
in a tight orbit around 2 And Ba, a late A or an early F star. 
The amplitude of the RV variation, Ki sini, would be equal 
to 10 and 20 km s”*^ for i = 90°and 30°, respectiv e ly. Ac - 
cording to the ephemeris provided by iRica Romerol (l2010l) . 
the separation of the components A and B will be 0.135”in 
2014 January, decreasing to 0.046”by 2018. In the near fu¬ 
ture, it will be thus impossible to get a spectrogram of 2 
And B without an overwhelming contribution from 2 And 
A. Consequently, measuring the RV of 2 And B may be very 
difficult. 

This paper has been typeset from a T]eX/ lAlj^X file prepared 
by the author. 
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